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ABSTRACT 



Using the finite-size effects the scaling dimensions and correlation functions 
of the main operators in continuous and lattice models of Id spinless Bose-gas 
with pairwise interaction of rather general form are obtained. The long-wave 
properties of these systems can be described by the Gaussian model with 
central charge c = 1. The disorder operators of the extended Gaussian model 
are found to correspond to some non-local operators in the XXZ Heisenberg 
antiferromagnet. Just the same approach is applicable to fermionic systems. 
Scaling dimensions of operators and correlation functions in the systems of 
interacting Fermi-particles are obtained. We present a universal treatment for 
Id systems of different kinds which is independent of the exact integrability and 
gives universal expressions for critical exponents through the thermodynamic 
characteristics of the system. 
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1 Introduction 



Recently the new effective method to investigate Id quantum systems using the finite-size 
effects in conformal field theory (CFT) was proposed in a number of papers [1-6]. In Id 
quantum systems there exist phase transition at zero temperature and corresponding 
long-distance conformal symmetry appears to be enough to determine the properties 
independent of the interaction at small scales. These properties are really of a great 
interest in the Id systems. Among them the correlation functions, i.e. expectation values 
of operator products in different points are especially interesting. The correlations at 
large distances decrease as a power of the distance (at zero temperature) the "critical 
exponents" being continuously dependent on the coupling constant [7-9]. 

In 2d systems the conformal invariance is known to constrain the spectrum of scaling 
dimensions [10]. It does really determine the critical exponents in Id quantum systems 
(i.e. in 2d models of the quantum field theory with one space and one time directions). 
The essential point is that it is possible to calculate central charge and scaling dimensions 
of the effective long-distance conformal theory using the so-called finite-size effects [1,2,11]. 

This paper is a review of applications of CFT methods to "real" physical Id systems 
which, as they appear, are neither conformal nor exactly integrable. We propose a unified 
treatment of physically relevant Id systems of different kinds. Among them are continu- 
ous and lattice models of Bose or Fermi multiparticle systems with interaction of a quite 
general type and Id lattice antiferromagnetic spin chains. We shall study the long-wave 
behaviour of the correlation functions in these models. The corresponding critical expo- 
nents are determined from the leading finite-size corrections to the low energy spectrum. 
A remarkable point is that these leading corrections are in a sense universal: they are 
completely determined by general thermodynamic characteristics of the systems. So our 
expressions for the critical exponents through certain susceptibilities have a wide range 
of applicability. They include all previously known examples as particular cases. 

This review is based mostly on our short publications [22,23,43]. Here we present the 
more detailed argumentation with new examples and discussions. 

To begin with, let us remind the ideology and main formulas of the finite-size approach 
in CFT [11]. Consider a CFT on an infinite strip of width L in space direction. Then each 
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primary conformal operator gives rise to an infinite "tower" of the eigenstates 10^) of 
the Hamiltonian with energies 



E*(k,k) = E?r+—(h+k+k) (i) 

and momenta 

O-rr 

Pt(k,k) = P++ T ( 8 + k-k). (2) 

(such that (vac\(p\(p k ^) ^ 0). Here h = A + A and s = A — A are scaling dimension and 
spin of the operator respectively, k, k > - are integers, E v ^ c is the ground state energy 
of the system in the "box" of length L. The momentum of the lowest state |0o,o) = 10) of 
the tower at L — > oo is denoted by P^. At last, the parameter v in (1) takes into account 
the difference in units of the space and temporal quantities, i.e. v is merely the velocity 
of sound excitations (group velocity at the Fermi surface). 

Thus, to find the spectrum of primary operators it is sufficient to calculate the energies 
of the lowest excitations Ef(0, 0) = of each tower up to the terms of order L^ 1 

E% - E v L ac = 2nvhL-\ (3) 

The states from 0-tower with k, k ^ correspond to the so-called descendants (or sec- 
ondary fields) of 0[1O]. They have the definite scaling dimensions h + k + k (but they 
are not primary fields!) and give a contribution to correlators too. However, we shall see 
that the leading term of asymptotics of the correlators is determined only by the primary 
fields. 

The long-wave asymptotics of two-point equal time correlator of the field has the 
form 

(0(z)0(O)> ~ cos(Pix)x- 2h , (4) 

the oscillating factor cos(P^a:) being determined by a gap into the momentum spec- 
trum. For simplicity, throughout this paper we consider only the equal time correlators. 
Their asymptotics is completely determined by the scaling dimensions h of the field 0. 
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The generalization to the space-time correlators is straightforward. In this case one also 
should determine the value of spin s. 

Indeed, the finite-size approach allows one to find out the central charge of the Virasoro 
algebra of the conformal symmetry using the equation similar to eqs.(l)-(3). Namely, the 
correction ~ L^ 1 to the ground state energy of the system is proportional to the central 
charge c [1,2]. For example, for the theory with periodic boundary conditions one obtains: 

7TC1) 

ET = e L - ™ (5) 

Here eo is the energy of the ground state per unit length in an infinite volume. The first 
term in (5) depends on regularization, but the second does not. In other words, the first 
(thermodynamic) term is determined by ultraviolet behaviour of the system; the leading 
corrections being given by the infrared properties. 

Usually it is convenient to calculate the central charge by means of another formula 
which is equivalent to (5). Due to the modular invariance the case of zero temperature 
(T = 0) and large L is equivalent to the case of the infinite volume (L = oo) and low 
temperature T. In the latter case (5) is rewritten as follows: 

vrcT 2 

f(T) = eo~, (6) 
where f(T) is the density of the free energy as a function of the temperature (at small 
T). 

In the paper [3-6] the equations (l)-(3),(5)-(6) were used to calculate the critical 
exponents in Id exactly solvable models. Namely, the low-lying excitation energies have 
been calculated within the Bethe-ansatz framework [12]. At present time there exist 
many papers devoted to this topic (see, for example, [1-6,13-20] and references therein) 
and we would not repeat here the finite-size calculations in integrable theories. Instead, we 
shall concentrate on the following important but less familiar aspects. First, the method 
turns out to be applicable to the models with general interactions and leads to results of 
the same completeness as in the exact solvable models. Second, the results which have 
been obtained earlier are in fact valid only for bosonic systems. The critical exponents 
in fermionic systems differ, generally speaking, from the bosonic ones and can also be 
obtained by means of the finite-size method. At last, the finite-size effects allow one to 



calculate the asymptotics of vacuum expectation values of some non-local operators in 
the Id systems. All these questions are the subject of the present paper. 

We shall consider the following classes of Id systems. The first one contains the 
continuous models of the spinless Fermi- or Bose-gas with interaction of a general form, 
the second quantized Hamiltonian being (the mass of particle is 1/2) : 

H = dx- d x i/j*(x)d x i(j(x) + \- 9 J J q dxdy ■ ip* (x)ip* (y)V (x - y)ip(x)i(j(y) . (7) 

Here L is the length of the system, V(x) is some even pairwise (repulsive) potential of 
a rather general form, g > is the coupling constant. The operators satisfy usual 

equal-time (anti) commutation relations (to indicate the statistics manifestly we sometimes 
use the notation ipB or ipp)- The number of particles in the system N is conserved; in 
the thermodynamic limit iV — > oo, L — > oo and p = N/L is the equilibrium density. 
Note that the potential V(x) = S(x) admits the exact solution [12]. In the first quantized 
framework the Hamiltonian (7) has the form 

# =-£^+s£n^-^)- (8) 

t=i ux i i<j 

Sometimes we use this language for convenience. 

The second case under consideration is the lattice version of (7)- (8): 

L 

H= - J2\-^( x+ l )^ x ) +^{x)%l)(x+ 1) -2^(x)ij(x)] 

x=l 
L 

+ 5 E ^(x)^(x) V(x - y) 4>\y)i>(y) (9) 

x<y 

Here we put the lattice spacing equal to 1, so x and L become dimensionless quantities 
now. If necessary, one can easily restore dependence of the lattice spacing in final answers. 

At last, we shall deal with XXZ Heisenberg antiferromagnet (or simply the spin chain) 
with the lattice Hamiltonian: 

1 L 

Hxxz = -r E^X+i + vlvl+i + COS7(7^ +1 ), < 7 < 7T. (10) 
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Here L is the number of sites of the lattice, are usual Pauli matrices, 7 is the anisotropy 
parameter. Isotropic (XXZ) Heisenberg antiferromagnet corresponds to 7 = 0. For the 
XXZ antiferromagnet model there exists an exact solution. Namely, all eigenstates of the 
Hamiltonian (10) can be constructed explicitly by the Bethe's method [12]. The simplest 
eigenstate with all spins looking at the same direction (i.e. full spin is equal to L/2) is the 
bare (unphysical) vacuum. The physical vacuum of the antiferromagnet has the minimal 
possible value of the full spin (0 or 1/2 depending on the parity of number of sites L) and 
corresponds to the filled bare vacuum, the number of reversed spins being an analog of 
the number of particles in the models (7)-(9). 

In fact, after the Jordan- Wigner transformation (that is, ip(x) = (1TJ=i a j)°'x > 4^{ x ) — 
c^"(n|=i CTj) (10) can be considered as a special case of (9). The simplest Hamiltonian 
(10) is chosen to make possible a comparison with the Bethe ansatz results. We can also 
consider more general spin 1/2 chains with the exchange interaction corresponding to 
arbitrary potentials in (9). 

All these models have a number of common physical properties. The most important 
property is the absence of a gap in the low-energy spectrum. In fact, such a general 
statement is not quite correct because one correctly can tune V(x) and p to create a gap 
in the spectrum of the lattice system (9). We are convinced, however, that the generic 
case with spinless Id systems is just the gapless spectrum with unique sound velocity^. 
Some comments on this point are contained in the Section 2. 

Another important property of models under consideration is that the central charge 
calculated from (4) in XXZ antiferromagnet is equal to 1 [3,14] and just the same holds 
for the model (7)-(9). So all these models are described by the universality class of a 
Gaussian model [21]. It is not very amusing since the lowest excitations are nothing but 
free phonons. 

Some comments on long-wave approximation are in order. The proper distance x c 
when individual particles have sense is equal to x c ~ L/N = 1/p. At large distances (x ^> 
x c ) there arises an effective theory describing the phonon excitations, the information 
1 As for the systems with internal degrees of freedom, they have in general a number of branches of 
gapless excitations with different sound velocities. This more complicated case is out of the scope of the 
present paper. 
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about interaction at small distances being contained in the only parameter v. Then the 
low-energy excitations corresponding to the linear dispersion law lead to effective scaling 
and conformal invariance (the 0(2)-invariance is trivially restored by the choice of units 
so that v — 1). Therefore, the deviations from the exact conformal theory can only be 
discovered at small L, so the series in L~ l are correct. 

Now we briefly describe the content of the paper. In the Section 2 we calculate the 
energies of low-lying excitations for the continuous Bose-gas model and obtain the spec- 
trum of scaling dimensions for this system. The more complicated case of lattice systems 
is considered in the Section 3. Here the universal expressions for the "compactification 
radius" in the Gaussian model through thermodynamic parameters of the system will be 
obtained. Then, using all these results, we write down the asymptotical series for corre- 
lation functions (Section 4). The most important examples (pair density correlator and 
one-particle density matrix) are discussed in detail. It turns out that a simple generaliza- 
tion of the finite-size method allows one to find the asymptotics of vacuum expectation 
values of some non-local operators as it is demonstrated in Section 5. Here the main 
examples are some special non-local operators in the spin chain. At last, the Section 6 is 
devoted to the systems of Fermi-particles. We show how the finite-size method may be 
used to calculate the exact critical exponents in this case. 

The Appendix is devoted to a non-trivial example of Id multi-particle systems with 
long-ranged interaction — the so-called Sutherland (or, Calogero) model [28]. First, we 
calculate the effective central charge from the low-temperature behaviour of the specific 
heat. Second, we compare the exact results for correlation functions in this model with 
the asymptotical series obtained from the CFT-approach. This model has never been 
discussed in the literature from the conformal point of view. It seems to us that the 
detailed exposition of this example could be useful. 

2 Scaling dimensions in the Bose-gas model 

Let us consider the Bose-gas described by the Hamiltonian (5). To obtain the scaling 
dimensions of different operators it is necessary to classify the low-lying excitations having 
zero energy at L = oo. We imply periodic boundary conditions, i.e. the particles live on 
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a circle of the length L. 

The typical excitation spectrum of the lei spectrum is shown schematically in Fig.l. 
One can learn that it is gapless with the energy tending to zero when the momentum is 
equal to 2irpm(m is integer). These properties hold for rather large class of both long-range 
and short-range potentials V(x). Indeed, the cases of small and large constants g in (7)-(8) 
have been studied in the literature. In [24], where the case of small g in Fermi systems was 
studied perturbatively, it was concluded that the correlators have a power-like behaviour 
and no gap emerges. The case of long-range V(x) and large g (strong repulsion) was 
considered in refs. [25,26,27]. In this case, the absence of the gap becomes quite obvious. 
Indeed, independently from V(x), the particles form a regular "dynamical lattice" (the 
so-called Wigner crystal) with sound-type low-energy excitations. From the other hand, 
there exist several exactly solvable models which possess a gapless spectrum. These are, 
for example, the model with delta-shape potential and that with the potential V(x) = x~ 2 
(the Sutherland model [28]). Evidently, the gap will not appear when performing small 
deformations of the Sutherland potential which do not entail a qualitative rearrangement 
of the ground state. 

We should investigate the spectrum when L is large but finite. Then the spectrum is 
quasi-discrete with the energy gap ~ L^ 1 (see Fig. 2). We are interested in the states with 
zero energy at L = oo. The states \<f>) corresponding to the primary conformal operators 
have the minimal energies in comparison with that of the neighbouring states. So each 
gapless branch of the spectrum "gives rise" to a primary operator. The higher excitations 
are the states of the conformal tower with /c>0or/c>0in(l). 

Certainly, we are able to calculate the energy E\ of the states |0) when the physical 
states can be constructed explicitly. This is the case for the exactly solvable models [12], 
in particular, for the Bose-gas with 5-shape interaction. But it turns out that in general 
case all that information in fact is not necessary and the energy of a "primary" state 
|0) may be obtained using simple thermodynamic reasoning. 

Here we shall demonstrate this for the system (7)-(8). Let us begin with the excitations 
which conserve the number of particles. The simplest suitable excitation is the creation 
of a phonon. The minimal possible momentum is ±27r/L (the dots A , A in Fig. 2). 
Evidently, the energy is equal to 2nv/L, and the gap in the momentum spectrum is 
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absent: P^ = 0. Then (1,2) imply h = ±s = 1. Let us denote (f)°± the corresponding 
primary operators. Then 

(0 O + (^)0° + (O)) = (0°»0°_(O)> ~ x- 2 ; (0° + (:r)0 o _(O)> = (0° (x)<(0)> = . (11) 

A more interesting excitation has the momentum P^ = 2np (A 1 in Fig.2). This excitation 

is produced by the "rotation" of the whole system with minimal angular momentum (i.e. 
it can be described by the first level of the "rotator"). In other words, this state can be 
produced from the ground state by transition to another inertial frame of reference. The 
periodic boundary conditions imply the quantization of the velocity of this new frame. 
More generally, one produces the family of the states \(fio,m), m is an integer (A 2 and so 
on in Fig.2). The momenta of these states are equal to lispm: P ,m — 27rpm, and their 
energies can be obtained by considering the uniform motion of the whole system, the 
momentum of each particle being equal to 2nm/L : 

*V=iV(Wi)'^^. (12) 

Li V 

Comparing (3) and (12) we find the dimensions of the operators o ,m : 

ho jTn = 2irpv~ 1 m 2 , m e Z; s = . (13) 
The above derivation may seem to be not quite correct since it neglects the quantum 

nature of the ground state. But more rigorous arguments lead to the same result. Indeed, 
let us turn to another inertial frame and use the first quantized language for convenience. 
Then the full (TV-particle) wave function is transformed as 

N 

^({rr,}) -> ^({a;,}) = exp(iq ]T x k ) ■ (14) 

k=i 

This function should be single- valued, so the momentum q is quantized: q = 2nm/L, 
m G Z. Now acting to \P by the Schrodinger operator (8) one obtains (12). 

The same result can be obtained also in the second quantization formalism. In this 
case one should work with the Hamiltonian (7) shifted by chemical potential term 
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H^H-fxN, (15) 
where N = Jq dxip^(x)ifj(x) is the operator of the number of particles, 



\i = lim 



Q E va 



L^oo ON 

p=const 

is the chemical potential. The Galilean boost operator is 



(16) 



K — / dxxip'(x)ip(x) (17) 
J o 

and the transition to another inertial frame can be described by 

H -> e~ iqk He iqk . 

We leave the calculation for the reader. 

The states \<j>± ), |0o,m) are excitations in the sector with the fixed number of particles. 
Now let us consider the excitations which change the number of particles. Now one should 
work with the modified Hamiltonian (15). In Bose case the adding of n particles results 
in the energy shift 

<o = \n\dEr/dN*) = f^eo/dp 2 ) = ^ ~ , (18) 

where eo again denotes the energy density. Here we used the well-known thermodynamic 
relation 

v 2 = 2p(d 2 e /dp 2 ) . (19) 
Now we can obtain the dimensions of the new family of primary operators 0„ 5 o : 

o 

VTL 

K,o = ^ s = - ( 2 °) 
Certainly it is possible to combine the above excitations, i.e. to add n particles and then 
to "rotate" the system at the m — th level. In this way we obtain the operators 4> n ,m with 
the dimensions 
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where we denote 



(21) 



R 2 = Birfw' 1 . (22) 

The dimensions /i ,i and /i 10 have been obtained first in [3] (in this paper the case of 
V(x) = S(x) was considered), where they were represented in terms of integrable systems 
theory. One can check that these dimensions coincide with ours, with the former being 
expressed through the sound velocity. It seems to us that the form of the answer given 
by (21)-(22) is more convenient, as it has sense for quite arbitrary potentials, and all the 
information about potential is contained in the unique parameter v (or R). In the next 
section we shall obtain a universal expression for R valid for lattice systems as well. 

In exactly solvable models one can calculate the central charge using the formula (2). 
Indeed, both the model with V(x) = S(x) in (7)-(8) and XXZ^spin chain (10) provide 
c = 1 [3]. In the Sutherland-Calogero model (V(x) = x~ 2 — [28]) the dependence f(T) 
at low temperature can be found using the results of [29], and (6) yields the answer 
c = 1 too (see the Appendix). A more general argument is the following. Since the 
scaling dimensions depend on a continuous parameter, the central charge should be not 
less than 1 [30] . Then one can always extract an "irreducible theory" having c = 1 from 
the initial theory [31]. On the other hand, the central charge greater than 1 implies the 
number of degrees of freedom greater than 1. However, it is impossible for the systems 
under consideration since we have assumed the gapless spectrum with the only branch 
of phonon excitations. The result c = 1 seems to be quite natural as just such a theory 
describes the phonon system in the long- wave limit. 

Note that the spectrum of scaling dimensions (21) also implies c = 1. Indeed, roughly 
speaking, almost all conformal theories with c = 1 are equivalent to the Gaussian models 
with the spectrum of dimensions like (21) which can be parameterized by the unique con- 
tinuous parameter (the whole classification of those theories would include the Gaussian 
models with corresponding orbifold line and three isolated models in moduli space of the 
conformal theories [24]). 
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More precisely, the Gaussian model is the 2d free massless scalar field theory with the 
action (here z = x + iy, z = x — iy) : 

Sr = 2tt / d2zdzL P d ~ zL P ( 23 ) 
which is evidently [/(l)-invariant. The field ip(z, z) takes values in the "circle" of radius 
R, i.e. one identifies ip and (p + 2nR. The spectrum of this theory is given by (21) [21,32], 
with the parameter R playing the role of the compactification radius in string theory. The 
operator 4 >0 +{4 >Q -) can be identified with the (anti)chiral [/(l)-current d z <p(dz<p), and 0„ i?n 
can be identified with the exponential of the free field: 

n>m ->: exp{ip(p + iptp} : , (24) 
where : . . . : denotes the proper normal ordering and the dimensions are 

(A, A) = ( P 2 /2,p 2 /2), /i = A + A, s = A-A, (25) 

(p, p) = (nR- 1 + mR/2, nR' 1 - mR/2); n, m G Z . (26) 

This theory is invariant under the dual transformation R — > 2/R, and R = \/2 is the 
self-dual point corresponding to the isotropic (XXX) antiferromagnet (7 = in (10)). 

3 Scaling dimensions in lattice models 

As we saw in the previous section, the continuous models have an exact symmetry which 
allows one to obtain the answer immediately, namely, they are invariant with respect 
to Galilean transformations. In the case of lattice models (9)-(10) the exact Galilean 
invariance is broken. Nevertheless, we shall see that it is possible to apply the finite-size 
method in this case too [43]. In this section the most general relations between critical 
exponents and thermodinamic parameters of the system will be obtained. 

Consider the system of spinless particles on Id lattice with the Hamiltonian (9) mod- 
ified according to (15). Then p = de^/dfi is the density of particles at zero temperature. 
We imply periodic boundary conditions in (9). It will be convenient to imagine the lattice 
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as being rolled up into a ring. Again, we shall assume that there is no gap in the energy 
spectrum of the system. In this case our arguments do not depend on the specific form 
of V(x). 

The operators 4>°± can be constructed in the same way as in the continuous case and 
their dimensions and correlators are given by the same formulas (11). 

Now we consider the analog of "rotating states" |0o,m)- By definition, |0o,m) is the 
eigenstate of H with momentum lispm and minimal energy. Let us represent it in the 
form (it is convenient to work in the second quantization formalism in this case) 

N N 
|0O,m) = ^ eX Pfen) J2 X j)^m(Xl, ■ ■ ■ ,X N ) J| ^(x k )\0) , (27) 

x 3 j + 1 fc=l 

where p m = 2nmL~ 1 , |0) is the bare vacuum, N is the (conserved) number of particles. 
When m = O|0 o ,o) = \vac) and \I> is the ground state coordinate wave function. We have 
extracted the exponential factor in (27) in order to make ^> m close to \l/o- This is really 
the case because p m ~ L^ 1 (note that in continuous systems ^ m = \l/o for all m). 

Acting to \(f>o, m ) by the Hamiltonian H (9) one can see that \l/ m should be the ground 
state coordinate wave function (in the iV-particle sector) of the following Hamiltonian 
H(p) &tp = p m : 

L 

H(p) = H + (l-cosp)^(^ t (x + l)^(a;)+^ t (x)V'(x + l) 

x=l 

L 

- isinj9^(^ t (x+ l)ip(x) - ip\x)ip(x+ 1) . (28) 

x=l 

Evidently, H(p) is produced from H = H(0) by the following transformation ip^(x) — > 
exp(ipx)ip^ (x) , ij}{x) — > exp(ipx)ip(x). Due to the gauge invariance it means that the 
uniform magnetic field B = p is applied to the system. In other words, H(p) describes 
the system (9) in the uniform magnetic field orthogonal to our ring lattice (Fig.3). Note 
that the magnetic flux through the ring is equal to 2nm (because p is quantized: p — pm). 

Passing from H to H(p) we can write for the ground state energy shift e m (up to L^ 1 ) 
the following general relation which is evident from the thermodynamic reasoning: 

1 <9 2 e 

SE 0>m = e m = -Lp 2 m (-g-^) p=0 = 2nvL~ 1 irrjv- 1 m 2 . (29) 
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Here eo is the minimal eigenvalue of H(p). When there is no Galilean invariance the 
current j = i(ip\x + l)ip(x) — ip^(x)ip(x + 1)) depends on momentum p according to a 
non-linear relation. So we have introduced the "current-momentum" susceptibility r\ in 
(29): 

= 8Peo = dj_ 
dp 2 dp 

Of course, the same expression for the energy shift (29) may be obtained by means of the 
perturbation theory, with exact value of the L~ x order term being given by the first two 
orders in p (see (28)). 

By comparing (3) and (29) we find the dimensions of the corresponding primary op- 
erators: 

h ,m = ■Kr]v~ 1 m 2 . (31) 

Another way to excite the system is to change the number of particles. Similarly to (18) 
we have in the Bose case 

8E^ = \ n *L-\^) = \n\-*L-\ (32) 

where 

d 2 e dp 

X= d7 = dH (33) 
is the usual susceptibility which is in a sense "dual" to rj. (Note that in general (19) is no 
longer true!) Comparing with (3) we obtain a new family of primary operators with the 
following dimensions: 

h n ,o = (47™ X rV . (34) 
Now let us note that the two susceptibilities rj, x are connected with the sound velocity 

v by a universal thermodynamic relation. The simplest way to see this is to consider 
the wave equation for long waves. Repeating standard arguments leading to the wave 
equation and taking into account the relation between current and momentum one has 
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d 2 u(x,t) d 2 u(x,t) 
X ^^~ = ^^x^ (35) 
(here u(x,t) is the displacement of medium and x should be considered as a continuous 

variable). Hence 

v 2 = VX~ L • (36) 
Therefore, combining the above excitations and making use of (36) we find the spectrum 
of dimensions in the form (21), where 

R 2 = Anvx ■ (37) 

This is the general expression for the compactification radius and (due to (21)) for the 
critical exponents. Note that in the case of the exactly solvable Hubbard model the critical 
exponents were expressed through the susceptibility x m the paper [44]. Let us consider 
some particular cases. 

If the continuum limit is performed in (9) the exact Galilean invariance is restored. In 
this case 

X=\ (38) 

V 

and (37) coincides with (22). 

In the case of the XXZ model (10) we have [14] (after re-interpretation of p, pi: p 
now is the total spin of the chain, ji is external magnetic field): 



v = , (39) 

7 



X = 9 , 7 s ■ • (40) 
27r(7r — 7) sin 7 



Hence 



R 2 = — (41) 

7T — 7 

which coincides with the Bethe ansatz results. The finite size corrections and scaling 
dimensions were studied in the framework of Bethe ansatz in papers [14,45]. 
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4 Asymptotics of the correlation functions 



In this section we shall obtain the correlation functions of the models (7)-(10) in the 
asymptotic series form. The most important correlators are density-density correlator 
H(x) = (p(x)p(0)), p(x) = ip*(x)ip(x) being the density operator, and the one-particle 
density matrix S(x) = {tp*(x)tp(0)). In the spin chain the similar objects are H(x) = 
(a^af) and S(x) = (<7+crf ), with a ± = a 1 ± ia 2 . 

The operators ip(x), p(x) have no definite conformal dimensions, since they do not 
behave properly under conformal transformations. Nevertheless, they are local operators 
and, as such, can be represented as linear combinations of the primary and descendant 
operators. For simplicity, we write this sum at the moment rather symbolically without 
accounting of the "descendant contributions" (their role will be discussed below). As we 
shall see, to obtain the leading asymptotics it is sufficient to take into account only the 
primary operator contributions. Then, the general expression for a local operator 0(x) 
has a form 

6(x) = £C,exp(iP*x)0(x), (42) 

where the sum runs over primary operators and are some numerical coefficients. A 
possible gap in the momentum spectrum leads to appearance of the oscillating factors 
exp(iP^x) in this formula. To find the correlation functions of the type (0(x)0(y)) it 
is sufficient to calculate pair correlators of primary operators <fi by using the well-known 
technique of conformal field theory [10,33]. The most convenient way to do this is to 
represent primary fields in the exponential form (24) and to average their products by 
performing the conventional Gaussian functional integrals with the action (23). 

The coefficients are non-zero if (vac\0(x)\<f>) ^ in the thermodynamic limit (here 
(vac\ is the physical vacuum; for brevity sometimes we merely use the brackets to denote 
the physical vacuum expectation value). For example, the results of the previous section 
yield the following series (for Bose statistics): 

oo 

i/>b(x) = c m ■ exp(2nipmx)(p ltm (x), (43) 

m=— oo 
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$b( x )= S C-exp(-27rimpx)0_i_ ro (x), 



(44) 



p( x ) = 0° ( x ) + 0° ( x ) + £ c ^ • exp(2irimpx)(j)o, m (x), 



(45) 



In the last formula o ,o denotes the identity operator. Using these series one can directly 
obtain asymptotical expansions of the correlation functions H(x) and S(x) as well as 
similar multi-point correlators. Note that the terms (0+(x)0+ (0)) = {4>°_ (x)4>°_ (0)) ~ x~ 2 
appear in the pp-correlator, with the corresponding exponent being equal to the canonical 
dimension of the operator p (independently of the coupling constant). The asymptotical 
series for H(x) always contains such term but in the most interesting examples this term 
is non-leading. 

Now let us discuss the descendant contributions. Generally speaking, the expression 
(42) should be supplemented by the descendant operators for each primary one satisfying 
the above selection rule. The conformal Ward identities imply that the descendant oper- 
ator 0(x) in a correlation function can be substituted by a number of differentiations of 
the same correlator in which the primary field (f>(x) stays instead of its descendant 0(x). 
Thus, the descendant contributions correspond to the terms with exponents exceeding 
the leading ones by a positive integer. For example, the first descendant of 0o,i gives a 
contribution to H(x) of the form x ~^ r2 ~ 1 sin(27rpx). So, accounting for the descendants 
carefully, we are write down the whole asymptotical series (x ^> p -1 ) : 



#(*)-p 2 = E^- 2 - fc +E^ 



-hm 2 R 2 



k=0 



m=l 



°° nk 
J2 A m,k cos(27rmpx + — )x~ 

k=0 1 



(46) 



m=0 



00 wk 
B m>k cos(2irmpx + — )x~ 

Lfc=0 1 



(47) 



Here A, B are some numerical coefficients. 

The formulas (46), (47) are in a complete agreement with some exact results. These 
are the calculation of the density matrix in impenetrable boson system (V(x) = S(x), 
g = oo in (7)-(8)) [34,35] and the exact result for the correlator H(x) in the Sutherland 
model at the special value of the coupling constant g — 4 [28]. These correlators do 
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really have the form of (46), (47) (Appendix). One can also compare these formulas 
with the answers of the papers [26,36], where the asymptotics of H(x) and S(x) was 
obtained by a direct calculation in case of the strongly repulsive particles with long-range 
interaction. In those papers the perturbation theory with respect to the small parameter 
g" 1 , where g is the coupling constant have been constructed. The result was that the 
parameter R should be determined by the same expression (22) provided that all orders 
of the perturbation theory are taken into account. However, R is small parameter when 
g is large, so only the "primary" contributions {i.e. the terms with k = in (46), (47)) 
were essential in refs. [26,36]. The rest terms are negligible on the background of slowly 
decreasing contributions ~ x -\rn 2 R 2 _ j us t the same reason does not allow one to obtain 
the term ~ x~ 2 in (46) by the method of ref. [26]. Meanwhile, there are models (e.g. the 
antiferromagnet with 7r/2 < 7 < ir) where this term becomes the leading one. 

The last check of (46), (47) is the calculation of non-equal-time correlators, both the 
charged and neutral ones, by a new exact method proposed in refs. [37]. To compare these 
results with ours one has to take into account spins of the operators and to consider the 
combinations x ± vt instead of x. If we suppose the symmetry with respect to exchange 
x + vt •<==>- x — vt, we reproduce the answers of papers [37,38]. 

Now we would like to discuss the meaning of the series (46), (47) with unknown coef- 
ficients. At the first glance quite arbitrary function can be represented in such a form. 
Nevertheless, the expressions (46), (47) contain an important information. 

First, it turns out that in all cases when the correlators can be obtained using other 
methods (see above) the coefficients A, B quickly decrease when k — > 00, so the series 
are well convergent. From the other hand, the above mentioned approaches as well as 
exact methods usually reproduce just several leading terms of (46), (47). Thus these 
expressions seem to be adequate to the problem. Second, the subseries corresponding 
to the descendants (the sums over k in (46), (47)) can be rolled up to smooth analytic 
function, but the sum over primary fields contains the fractional degrees of x : 

H(x) - P 2 = ^+t f m {x) ■ x~ m2R2 l\ (48) 

X m=l 
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S{x)=Y,im{x)-x- 2 ' R2 - m2R2 l\ (49) 

m=0 

When R 2 is irrational, the number of different types of ramification is infinite (x a and 
or give the same ramification if a and (3 have the same fractional parts). It apparently 
implies the essential singularity. On the contrary, if R 2 is rational the correlators (48), (49) 
contain the finite number of the ramification types. This case corresponds to the rational 
conformal theories [39], and demonstrate the finiteness of the number of fields which are 
primary with respect to some large chiral algebra [40]. 

The question emerges whether the models described by the rational conformal theories 
in the long-wave limit have any peculiar properties. Our conjecture is that the ground 
state wave function in such models is analytic function of every variable on the finite cov- 
ering of the complex plane (may be with puncture at oo). Using the results of refs. [28,29] 
one can easily check that this is the case for the Sutherland modelQ. 

Now we would like to discuss briefly some corrections to the expressions (46), (47) 
originated from the deviations of the critical point Hamiltonian from the fixed point one. 
That is they can differ by the irrelevant operators which have their conformal dimensions 
hi rre i greater than 2 (and, therefore, do not effect the critical behaviour). Then they 
induce contributions to (46), (47) decreasing with the distance faster than the leading 
terms being suppressed by x irrel+2 . Though the structure of series in this case is more 
complicated their global structure is the same in main features. In particular, the notion 
of rationality (the finiteness of the ramification types) remains unchanged. Certainly, in 
special cases the corrections merely ought not to appears as it is the case for impenetrable 
Bose-gas. 

More interesting point is the possible logarithmic corrections to leading asymptotics. 
The reason of this is the appearance of the additional marginal operators (i.e. operators 
with dimension (1,1)) which might be given rise in critical Hamiltonian at some special 
values of coupling constants (and, therefore, R). As a rule, the logarithmic corrections 
emerge when the symmetry of the system increases [32]. For example, in the SU(2)- 
symmetrical case of isotropic antiferromagnet in zero magnetic field the leading term of 
2 We thank M.Olshanetsky and A.Perelomov for the discussion on this point. 
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the correlator (a x ai) given by (46) should be multiplied by (Inx) % [41]|. In less symmetrical 
models the multiplicative corrections to leading terms are absent, but they can appear in 
sub-leading terms. 

To conclude this section two remarks are in order. 

One can easily obtain the multi-point correlators using the expressions (43)-(45), as 
primary fields can be realized as exponentials of free scalar field and can be averaged 
directly with the Gaussian action (c = 1!). The final expression has a type 

(II P( x i)) ~ X! ° m i exp(i27rp^mj) \x{ — Xj\~^ mimjR2 + non-leading terms (50) 

» {mi} i<j 

and should be compared with the one which can be obtained by the method of [36]. Then 
the primary field contributions are the same in both cases (descendant terms, as above 
mentioned, can not be taken into account by this method). 

The second remark concerns the coefficients in the sums (43)- (47). Generally speaking, 
the arguments of trigonometric functions in (46), (47) should be added by some constant 
(or slowly varying) phase (if c m ^ c_ m , c* m ^ —c*_ m , c' m ^ c'__ m and so on). Nevertheless 
sometimes this phase is absent. For example, this is the case for (46) since c' m = c'_ m due 

to the evident symmetry with respect to the exchange q < > —q in the wave function \I/ 

of the other inertial frame. Here it would be appropriate to remark that the exponential 
operators </>o,m + 0o,-m from (16) form closed operator algebra which corresponds to global 
0(2)-invariance instead of £7(l)-invariance as it is in (23). 

Thus, we have demonstrated that the long wave correlation properties of bosonic 
spinless systems can be described by the Gaussian model with proper "compactification 
radius" R. The parameter R is defined by the general formula (37) which gives (22) and 
(41) in the cases of continuous bose-gas and Heisenberg antiferromagnet respectively. 

5 Vacuum expectation values of non-local operators 

In the main body of this section we shall consider the Heisenberg spin chain (10). Let us 
consider the following non-local operators in this model: 

3 We are indebted to V.E.Korepin who drew our attention to this paper. 
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y 

S(x,y) = J2 a j = ex P{^(^,y)}, (51) 

j=x 

q(x, y) being the operator of the number of reversed spins at the sites on the way from x 
to y, and 

T(x, y) = P(x, x + l)P(x + l,x + 2)...P(y-l, y)P(y, x), (52) 

where P(x,y) = 1/2(1 + a x a y ) is the operator of permutation of the spins at the sites 
x and y. The operator T(x,y) is simply the operator of cyclic permutation of the sites 
along the segment [x, y] : x — > x + 1, x + 1 — > x + 2, . . . , y — 1 — > y, y — > x. 

Our aim in this section will be to find the asymptotics of the following correlators 
(vac\S(x, y)\vac), (vac\T(x,y)\vac) , {vac\S(x, y)T(x, y)\vac) at \x — y\ ^> 1. Such non- 
local correlators are quite interesting objects from the point of view of the quantum 
inverse scattering method [42]. Another application of our results is the calculation of 
correlation functions in the system of fermions with spin 1/2 [27] where one should know 
the asymptotics of the above non-local correlators. Due to translational invariance these 
correlators depend only on \x — y\. 

We shall show that these correlators have a power-law asymptotics. The corresponding 
exponents can be calculated using a natural generalization of the finite-size method. The 
idea is the following. As it was shown in [21], the (extended) Gaussian model has a 
non-local sector in which one should choose the numbers n,m in (21) to be half-integers. 
The corresponding operators, e.g. o i, are so-called disorder operators. They can be 
expressed through the free field ip (see (23)) only in a non-local way. Acting at the 
ground state of the extended Gaussian model by such an operator, say 0o,±> one obtains a 
state in the sector with antiperiodic boundary conditions (b.c). Therefore, the correlators 
of even number of disorder operators have sense in the sector with periodic b.c. A couple 
of disorder operators at the points x, y in the spin chain looks like a non-local operator 
acting in the segment [x,y]. Certainly, that is what allows one to find the asymptotics of 
the vacuum expectation values of non-local operators by CFT. 

Let us consider an extended Hilbert space containing all states of the spin chain with 
different number of sites and different b.c. simultaneously. We introduce the following 
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operators acting in the extended space: the operator a/^{x){(3 = ±1/2) creating a new site 
of spin (3 between the sites x and x + 1 of the initial chain; the operator bp(x) annihilating 
the site x (of spin (3x) in the case of (3x — (3 and giving zero if f3x = —(3. Thus, ap(x) 
acts from the sector of the extended Hilbert space corresponding to the chain with L sites 
to the sector corresponding to L + 1 sites. Quite similarly, bp(x) changes the number of 
sites from L to L — 1. The operator b±i can be identified with non-local operator ± i 
in extended Gaussian model above. We also introduce the operator 

s(x) = n oj, (53) 

j=x 

which connects the sectors corresponding to periodic and antiperiodic b.c. 

Indeed, let us consider the ground state wave function ^(xi, . . . ,%) of the antifer- 
romagnet. Here Xk are the coordinates (integer numbers) of reversed spins. For periodic 
b.c. one has *&(xi) = ^(xi + L) with other Xk being fixed (for brevity we write down 
explicitly only one variable). Let us denote \P = s(x)^> . In order to find ^f(xi + L) we 
shall move X\ and keep other variables fixed. Suppose we have k reversed spins to the 
right from the site x. Then one can write for x\ < x : 

= s(x)f (xi) = (-l)^(x!). (54) 
For x\ > x one has the following evident equalities: 



f (xi+L) = s(x)^( Xl +L) = (-l) k+1 y( Xl +L) = (-l) k+1 ^( Xl ) = -s(x)^( Xl ) = 

(55) 

So ^ belongs to the sector with antiperiodic b.c. 
Now we can write 



T (^> y) =J2 a U x ) b p( x )' ( 56 ) 

5(x,y) = s(x)s(y). (57) 

Let us denote by \L, +)(\L, — )) the ground state of the spin chain with L sites and periodic 
(antiperiodic) b.c. (in the thermodynamic limit these states coincide with \vac), but we 
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are interested in finite size corrections ~ L~ v ). Working with extended Hilbert space we 
can consider the ground states in the sectors with another L and antiperiodic b.c. as 
excited states over \vac) = \L,+). Roughly speaking, a^(x) creates an excitation with 
spin 1/2 ("one half of magnon"). If L is even then only even number of such excitations 
can exist, and if L is odd then there necessarily exist odd number of such "half-magnons" . 
Clearly, in the latter case the ground state is two-fold degenerate. 
One can easily check that 

(L - 1, +\b s (x)\L, +1) + 0, (L + 1, +\a\(x)\L, +) + 0, (L, -\s(x)\L, +) + 0, (58) 

in the thermodynamic limit, i.e. the states \L ± 1,±) satisfy to the selection rule 
of Section 3. This means that in order to find the asymptotics at \x — y\ ^> 1 of 
(vac\T(x,y)\vac) = 2(a^i(x)6 + i(y)) (obviously, (a\ib,i(y)) = (a^i (x)b_i (y))) and 
(vac\S(x,y)\vac) = (vac\s(x)s(y)\vac) we can use the formula (4) with the scaling di- 
mensions being determined from the relation (3). It is sufficient to calculate the energy 
shifts of the states \L ± 1, +) and \L, — ) up to the first order in L^ 1 and find their mo- 
menta. This can be done using the well-known exact solution (Bethe-ansatz) [12,14]. 
Certainly, one should modify the Hamiltonian as follows: Hxxz — > Hxxz — toL where eo 
is the energy per one site. We omit these lengthy calculations here. 

The results are as follows. The spectrum of scaling dimensions is given by the formula 
(21) with half-integer n,m. The parameter R is determined according to (41). The 
correlation functions look like 

< vac\S(x, y)\vac >=< vac\ exp{i%q(x, y)}\vac >~ cos[7r(a; — y)/2] ■ \x — y\~ x , (59) 

< vac\T(x,y)\vac>~ \x - y\~^, (60) 
< vac\T(x, y) exp{i7rg(x, y)}\vac >~ cos[7r(a; — y)/2] ■ \x — y\~ x ~^. (61) 

where 
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// = R~ 2 /2, A = R 2 /8. (62) 

Note that in the case of isotropic (XXX) spin chain (7 = 0) the operator s(x) has scaling 
dimension 1/8 and coincides with the spin field [32]. 

The correlators (59)- (61) appear in the calculation of correlation functions in the 
multiparticle spin 1/2 systems as auxiliary objects. In particular, the result (61) was used 
in [46] for finding the long-wave asymptotics of the density matrix in the system of spin 
1/2 Fermi particles with strong interaction. 

To conclude this section let us make a remark on the corresponding non-local operators 
in the Bose-gas model (7)-(8). An analog of S(x,y) is defined by the same formula 

S(x,y) = exp[iirq(x,y)], (63) 

where q(x, y) is now the operator of number of particles in the segment [x, y]. The leading 
term of (S(x,y))(\x — y\ p _1 ) can be obtained using the above technique: 

< vac\S(x, y)\vac >~ cos(np\x — y\) • \x — y\~ R2 ^ 8 . (64) 
It is convenient to introduce an analog of the operator s(x) (53): 

s(x) = exp[inq(x, L)], (65) 

which will be useful in the next section. The meaning of the operator T(x,y) in the 
continuous case is obscure. 

6 Critical exponents in fermionic systems 

So far, dealing with the model (7)- (8), we implied the case of Bose-statistics. It is known 
that the critical exponent of the fermionic field correlator S F (x) = (ip F (x)ip F (0)) differs 
from that of the bosonic one Sb(x) = {i/j b (x)^b(0)) [26]. In this section we shall obtain 
the answer in the framework of the finite-size method. 

There are two different approaches to this problem. The first one makes use of the 
fact that the operator s(x) (65) produces the Jordan- Wigner transformation from bosonic 
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field operators to the fermionic ones: 

ip F (x) = ip B (x)s(x), 

r F (x) = s*{x)r B {x). (66) 
Following the line of the previous section we obtain the result: 

S F (x) ~ cos(vrpx) • x- 2,R2 - R2/ \ (67) 

It is also easy to write down the whole asymptotical series for Sp(x) and see that it is in 
a good agreement with that of the paper [36] obtained by completely different methods. 

Less formal arguments are the following. Let if>(xi, . . . , xjy) be the ground state wave 
function of the system (7) or (8) (in the first quantization language). It should be an 
eigenfunction of the full momentum P with an eigenvalue P. So we can write 

exp(iPa) • *&(xi, . . . , xn) = ^{x\ + a, . . . , xpj + a) = exp(iPa) • \l/(xi, . . . , x^), (68) 

where a = L/N. On the other hand, the cyclic permutation Xi — > results in the sign 
factor: 

*(x 2 , x 3 ,..., x N , Xl ) = (-lf-^ixu ...,x N ). (69) 

Setting xi = x, X2 = x + a, X3 = x + 2a, . . . , xn = x + (N — l)a and comparing these 
formulas^ we obtain the following selection rule for the full momentum: 

exp(tPL/N) = (-lf-\ (70) 
or, writing P = 27rpm, as before, 

(-l) 2m = (-l)^- 1 . (71) 

Here we use the fact that the ground state wave function ^(x, x + a, x + 2a, . . .) is non-zero everywhere 
(as it have no nodes at all). 
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Thus, when N is even, the ground state is two-fold degenerate (P = ±7rp) . In the case of 
free fermions this fact is obvious. 

Clearly, the neutral excitations with zero energy ("rotations" of the system; see Sect. 2) 
should have momenta P + 2nmp(m G Z) where P is the ground state momentum. Thus, 
the dimensions of o ,m are given by the same formula (13). It is easy to see that this 
is also the case when we add even number of particles to the ground state. To avoid 
misunderstanding it should be stressed that the dimensions in bosonic and fermionic 
cases are generally not equal to each other, the point is that the sound velocity may be 
different in the two cases. 

But the addition of an old number of particles n changes the parity of N, so, according 
to (71), we have for n — 1 : 

5E[ fl = 2ttvL~ 1 (R' 2 + R 2 /1Q) = 2-KvL~ 1 h l i_ ) (72) 

where we have formally used the notation from (21). 

Combining the different excitations in accordance with (71), we obtain the general form 
of scaling dimensions in the fermionic systems. The spectrum is given by the formula (21) 
with the following conditions on n and m : 

a) when n G 2Z , m G Z; 

b) when n G 2Z + 1 , m G Z + 1/2, 

i.e. m — m'-\ — , m! is integer. (73) 

7 Concluding remarks 

In this paper we have demonstrated how one can find the long-wave asymptotics of various 
correlation functions in one-dimensional field-theoretical models using the finite-size tech- 
nique in CFT. The obtained results are quite general. They are valid both for integrable 
and non-integrable models with interaction of a general form. 

In conclusion, we would like to say a few words about the meaning of our results. We 
have shown that the long-wave correlation properties of the wide class of one-dimensional 



26 



spinless systems are described by the Gaussian model, the simplest conformal field theory. 
The Gaussian model still contains a free parameter — "compactification radius" R. 

An important point is that the value of R is determined by the nature of the short-range 
interaction (of the order of the mean distance between the particles x c ). So the short- 
range properties of the system are included into the effective long-range CFT through the 
only one parameter. 

One of the main results of this paper is the expression for R in concrete systems through 
their thermodynamic quantities. Clearly, it immediately gives us all the critical exponents 
as functions of thermodynamic parameters of the system. At last, for a given system it is 
much easier to solve the thermodynamics than to calculate correlation functions directly. 

We are grateful to V.E.Korepin, A.Yu.Morozov and A.A.Ovchinnikov for fruitful dis- 
cussions. 



8 Appendix 

Let us consider as an example the system of the type (8) with 



V(x) = x- 2 (Al) 



or 



V{x) = ^J^TT, (A2) 



sm 2 (27rx/L) 

(in the thermodynamic limit (Al) is equivalent to (A2)). The model (Al) is called the 
Calogero model, (A2) — the Sutherland model [28] . We call its thermodynamic limit the 
Sutherland-Calogero (SC) model. 

In order to find the effective central charge we should investigate the behaviour of 
the free energy f(T) at low temperature T — > 0. The thermodynamics of SC model was 
studied in [29]. 

First, we have the general thermodynamic relation 



f(T)=w-P, 
27 



(A3) 



where P is the pressure. In [29] it was shown that the thermodynamics of SC-model is 
described by the system of equations 



^wr^ra- « A4 » 



^^r* 1 ^' <a5) 

where (3 = T™ 1 , a = (1 + y/1 + 2g)/2 (g is the coupling constant in (8)), 



<j>(x) =/3fi + x + (1 - a) ln(l + e~ x ) . (A6) 

For finding power-like corrections in j3~ l at (5 — > oo we can put x = — oo and expand 
sqrt(f> in a power series. Let 

/oo r , n 

dx . - ■ . 2 x + (1 - a) ln(l + e"*) , (A7) 
-oo [c + ct J L J 

then we rewrite (A4) up to f3~ 2 as follows: 

i al\ Q.I2 , . oX 

anp = mu 2 H 5 3- , (A8) 

2/^3 2/3 2 /J 

3 3 Sot/jiIi 3al2 

-anp = run* + + ——^ . (A9) 

The computation of the integrals gives 

h = , h = y , (A10) 



then 



The sound velocity is [28,29] 

v = 2nap , (A12) 
so, comparing with (6), we conclude that c = 1. 
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The exact density-density correlation function H(x) is known in this model at g — 4 
[28]: 



1 r 2x 

H(x) = F 2 (2x) - -F'(2x) / dzF(z) + 1 , (A13) 

2 jo 

where 

„, sin7ra; , k 

F(x) = (A14) 

(we put p = 1). In this case a = 2 and i? 2 = 2. 

Expanding (A13) in a power series of x~ k we obtain exactly the series (46) which 
follows from the CFT. 
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Figure Captions 

Fig.l The excitation spectrum of Id spinless Bose-gas (schematically) 
Fig. 2 The spectrum of the same system in a finite volume 
Fig. 3 The lattice system in the uniform magnetic field 
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